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Abstract 

The Foulkes conjecture states that the multiplicities in the plethysm 
Sym“(Sym^y) are at most as large as the multiplicities in the plethysm 
Sym*’(Sym“y) for all a < b. This conjecture has been known to be true 
for a < 4. The main result of this paper is its verification for a = 5. This 
is achieved by performing a combinatorial calculation on a computer and 
using a propagation theorem of Tom McKay from 2008. 

Moreover, we obtain a complete representation theoretic decomposi¬ 
tion of the vanishing ideal of the 5th Chow variety in degree 5, we show 
that there are no degree 5 equations for the 6th Chow variety, and we also 
find some representation theoretic degree 6 equations for the 6th Chow 
variety. 
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1 Introduction 

The Foulkes conjecture |Fou50) states an inequality of certain representation 
theoretic multiplicities. For a finite dimensional vector space V, let Sym^U 
denote its 6th symmetric power and let Sym“Sym^U denote the ath symmet¬ 
ric power of Sym^U. We will always assume that dimU is large enough, i.e., 
dimU > max{a, 6}. As a GL(U) representation, Sym“Sym^U has an intriguing 
structure with interesting long-standing open questions like the following one. 

Conjecture 1 (Foulkes conjecture). Let a,b € N>o with a < b. Then for 
every partition A the multiplicity of the irreducible GL(K) representation {A} 
in the plethysm Sym“Sym^K is at most as large as the multiplicity of {A} in 

Sym^Sym“y. 
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The inequality a < b is important: We know that Sym“Sym^y contains 
irreducible GL(y) representations {A} for which the Young diagram of A has up 
to a rows, but Sym^Sym“Y contains irreducible GL(Y) representations {A} for 
which the Young diagram of A has up to b rows. 

Brion [Bri93) showed that Conjecture [T] is true in the cases where b is large 
enough with respect to a. Conjecture [1] is true if we only consider partitions A 
with at most 2 rows |Her54] , a phenomenon called Hermite reciprocity. Manivel 
showed that Conjecture [T] is also true for all partitions A with very long first 
rows [Man98[ Thm. 4.3.1]. In this case we do not only have an inequality of the 
multiplicities, but equality. 

Conjecture [T] is known to be true for all a < 4: For a < 2 see the explicit 
formulas in |Thr42) : for the case a < 3 see |DS00) . The proof for a = 4 uses a 
canonical map 4/4 4 and a propagation theorem by McKay. This study is related 
to a stronger conjecture, sometimes called the Foulkes-Howe conjecture, which 
at least dates back to Hadamard |Had97| . We will explain this in more detail 
now. Clearly Conjecture [1] is equivalent to saying that there exists a GL(I4) 
equivariant inclusion map 

Sym“Sym^I4 ^ Sym^Sym“I/. (2) 

A natural candidate for the map in ([2]) was the following map d'a 

Sym“Sym^K —Sym^Sym“K 

r\ 

L Q 

{g)“{g)''K-^^ {g)^{g)“l/ 


where t denotes the canonical embedding of symmetric tensors in the space of 
all tensors, g is the canonical projection, and r is the canonical isomorphism 
given by reordering tensor factors. 

Hadamard conjectured [Had97) that '^a,b is injective for all a < b. Howe 
|How87[ p. 93] wrote that this “is reasonable to expect”. This conjecture is 
known as the Foulkes-Howe conjecture. However, Muller and Neunhoffer |MN05j 
showed that 'I's 5 has a nontrivial kernel. As is in our case, the results in |MN05j 
were also based on computer computations. 

McKay [McK08] contributed the following important propagation theorem. 

Theorem 3. If 4'^ is injective, then 4>a,c is injective for all c>b. 

Theorem[3]allows us to verify Conjecture[T]in infinitely many cases while only 
doing a finite calculation: |MN05j calculate that 4 / 4 x 4 is injective, so Theorem[3] 
implies the following corollary. 

Corollary 4. ConjectureUi is true for a = 4. 

Our main result is the following theorem. 

Theorem 5. Conjecture{li is true for a = 5. 
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We obtain this result as a consequence of the first part of the following 
theorem. 

Theorem 6. (a) 4'5_6 injective. 

(b) As a G\-{y)-representation the kernel o/^'s^s decomposes as 

ker vI/ 5.5 = {(14, 7, 2, 2)} © {(13, 7, 2, 2,1)} © {(12, 7, 3, 2,1)} © {(12,6, 3, 2, 2)} 
© {(12,5,4,3,1)}© {(11,5,4,4,1)}© {(10,8,4,2,1)}© {(9,7,6,3)}. 

(c) Let S denote the set of partitions A that satisfy 

• A has 36 boxes, 

• A has 6 rows (i.e., Ag ^ 0} and 

• A appears with multiplicity exactly 1 in Sym®Sym®14. 

The following types appear in ker'I'6,6 (of course with multiplicity 1), and 
no other partition from S appears in ker'I'6,6' 

(20,7,6,1,1,1), (9,8,8,7,3,1), (16, 7,5,3,3,2), (15,9,4,3,3,2), 

(11,6,6,6,4,3),(11,9,7,3,3,3). 

Although 'I' 5,5 is not injective, if a = &, then Sym“Sym^y = Sym^Sym“I4. 
Therefore Theorem and Theorem [ 6 }a) together imply Theorem 

Remark 7. Regarding Sym^R as polynomials over V, let ‘ift, := 

G\-v{X 1 X 2 • ■ • Xh) C Sym^R denote the GLy orbit closure of the monomial 
X 1 X 2 ■ ■ ■ Xb- We call this orbit closure the 6 th Chow variety. The kernel of 'i>a,b 
is the homogeneous degree a part of the vanishing ideal of '^b, as was shown by 
Hadamard (see e.g. |Lanll[ Section 8 . 6 ]). This interpretation is what greatly 
inspires the algorithm described in this paper, which is used to obtain the results 
in Theorem [51 

Moreover, we remark that ^b is an important subvariety of the determinant 
orbit closure in geometric complexity theory |Lanll] , |Kuml2] , [BHI15] : an ap¬ 
proach towards resolving major open questions in theoretical computer science 
via geometry and representation theory. 

Theorem | 6 ] can be interpreted as follows: There are no equations for “^6 
up to degree 5, the vanishing ideal of in degree 5 is identified as a GLy- 
representation, and some parts of the vanishing ideal of in degree 6 are 
determined. ■ 

For more information on the history of the Foulkes conjecture we refer the 
interested reader to [Lanl51 Section 7.1]. 

We note that the result in Theorem [ 6 }b) uncovers the precise structure of 
ker'I' 5 , 5 , which is slightly different from the result “with high probability” in 
[Ti^’ (6.1.4)]. 
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Organization of the paper 

As stated above, Theorem [S] follows from McKay’s theorem and Theorem [SI 
The rest of the paper is devoted to the the algorithm behind the computations 
in Theorem HI In Section |2| we use representation theoretic considerations to 
reduce the problem of computing ker'f'o b to computing the dimension of the 
kernel of the restriction of 4'a,b to the subspace of Sym“Sym^K of highest weight 
vectors of type A for all A. In Section |3| we give a detailed description and a 
justification of the algorithm to compute this dimension. An overview is given in 
Section l3.ll While the algorithm includes various mathematical optimizations, 
the problem is still computationally too big to be feasible. In Section 0] we 
describe various techniques we used to speed up the algorithm and make the 
calculation feasible. In Section |5| we give remarks on the implementation and 
running time of the algorithm. 

Acknowledgments 
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valuable discussions. This work was initiated at the 2012 Mathematics Research 
Communities Geometry and Representation Theory Related to Geometric Com¬ 
plexity and Other Variants of P v. NP in Snowbird, Utah. We thank Texas 
A&M University for its hospitality for a one week collaboration effort funded 
by an AMS MRC Collaboration grant. 

2 Reduction to highest weight vectors 

Our goal is to compute the kernel of the GL(U)-morphism and decompose 
it as a GL(U)-representation. Using Schur’s lemma we can treat each GL(U)- 
isomorphism type A independently and calculate the kernel of the restriction 
the A-isotypic component of Sym“Sym^U: 

kervl/,,b = 0ker'k^;r, 

A 

where the sum is over all partitions A with |A| := = ab. The following 

definition simplifies our computation. 

Definition 8. Let U C GL(U) denote the group of upper triangular matrices 
with Is on the main diagonal. We denote by diag(ai,..., Odim v) the diagonal 
matrix with complex numbers oi,... jOdimV on the main diagonal. Let W be 
a polynomial GL(U)-representation. Let A be a partition. A vector f € W that 
is invariant under the action of U and that satisfies 

diag(ai,..., Odim y) • / = • • • adto”y / 

is called a highest weight vector of type A. ■ 
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Let W := Sym“Sym^F and let W := Sym^Sym“y. Let HWVa(VL) denote 
the vector space of highest weight vectors of type A in W. Since every irreducible 
GL(y)-representation has a unique highest weight vector (up to scale), we only 
have to check the kernel of the restriction '1'^^ : HWVA(kL) —S' HWVA(kL') of 
to HWVA(fL): 

ker«'^j° = linspan{g/ | g S GL(F),/ e ker^-^ J. 

So kerd>^“),®° is a A-isotypic GL(1/) representation that contains the irreducible 
representation of type A with a multiplicity of dimker'k^ j,. We write: 

multA(ker= dimker'I'^j. 

Hence we are left with determining dimker 

Let p := dimHWVA(kL) and p' := dim HVVVa(VL 0- cases we have 

p < p' (otherwise this would disprove Conjecture [T]) . To compute dimker j, 
we first determine a basis {/i}i<i<p of HWVa(W) and a basis {fi}i<i<p' of 
HWVa(W'). Using these bases we can write down the p' x p transformation 
matrix of and determine its kernel with linear algebra. 

3 The algorithm 

3.1 Overview 

In our cases for every A the natural numbers p and p' can be readily com¬ 
puted using the computer programs Schui0 or Li^. We then construct a basis 
of HWVa(IU) by choosing more and more “random” elements in 
HWVa(IU) and by checking their linear independence until we have p many lin¬ 
early independent ones. In the same way we construct a basis {/{,..., /p/} of 
HWVa(IU'), but this time we construct it together with a set of p' distinct points 
{v[, ... ,Vp,},Vj G Sym“U, such that the p' x p' matrix of function 

evaluations has full rank. This means that every element /' of HWVa(IU') is 
uniquely defined by its evaluation list ..., f'{Vp,)). Note that here the 

fi G Sym^Sym“U are regarded as functions on Sym“U as opposed to vectors 
from the tensor product space. 

Then we apply I'a.b to /i,..., /p and determine the rank of the p' xp matrix 
given by the function evaluations ■ The kernel of is given 

precisely by the kernel of this matrix. In other words a linear combination 
Yh=i cti'^aAfi) vanishes identically iff J2i=i = 0 for all j. 

^online available at http://sourceforge.net/projects/schur 
^online available at http://www“inath.univ-poitiers.fr/-'inaavl/LiE 
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3.2 Numerical considerations 

Even though at a first glance the method described in Section IXT] might look like 
it works only “with high probability” or “numerically”, it gives symbolic results. 
The key property that we use is that the p' x p' matrix has full rank 

by construction, so that a function in HWV;,(VE') is uniquely determined by its 
evaluations at the points 

We will construct the functions fi, f- over the integers and we will also do 
so for the points u'. Therefore all evaluation results are integers. Integer calcu¬ 
lations can be handled without any numerical errors by software with bignum 
arithmetic support. 

3.3 Constructing the highest weight vector bases 

In this section we describe how to obtain the functions /i, 1 < i < p and //, 
I < * < p'. Since their number of monomials is very large, we will not compute 
the monomials, but we will provide a method to evaluate the functions without 
writing them down. This is the same phenomenon as for the classical n x n 
determinant, which has n\ summands but which can be computed in polynomial 
time. 

3.3.1 A basis for the tensor power 

A Young diagram is a left- and top-aligned array of boxes. We identify partitions 
with Young diagrams by interpreting the entries in the partition as row lengths. 
For example, the Young diagram corresponding to the partition (4, 2) is 


A Young tableau is a Young diagram in which we write a natural number in 
each box. Different boxes can have the same entry. A Young tableau is called 
semistandard if the numbers increase from top to bottom in each column and 
never decrease from left to right in each row. An example of a semistandard 
Young tableau is the following: 


1 

1 

3|3 

2 

2 



If a Young tableau T with ab boxes has exactly b Is, exactly b 2s, ..., and 
exactly b as, then we say that T has rectangular content a x b. The example 
above has rectangular content 3x2. A special case are the Young tableaux 
with d boxes and rectangular content dx 1: Each number from 1,... ,d appears 
exactly once. A semistandard Young tableau with rectangular content d x 1 is 
called a standard tableau. 

We want to construct a basis of HWVA(Sym“Sym*'E). There is a nice basis 
of HWVA({^“Sym*'E) that is obtained from a nice basis of HWVa({^“^E) by a 
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process we call inner symmetrization. We describe this basis of HWVa({^“^ 1^) 
now. Let {ei}i<i<dimV be the standard basis of V. We start with a canonical 
construction of a highest weight vector of type A in For each column 

of length fc in A we take the wedge product ei A 62 A ■ • • A and tensor them 
together from left to right to obtain the tensor (x G For example, for 

the partition (4, 2) we get 

C(4.2) := (ei A 62) (81 (ei A 62) (g) ei (8> ei (10) 

and for the partition (7, 3, 2) we get 

C(7.3.2) := (ei A 62 A 63) (g) (ei A 62 A 63) 0 (ei A 62) 0 ei 0 ei 0 ei 0 ei. 

The vector (^\ is a highest weight vector of type A in which is easy to 

verify (recall Definition |S|). The group &ab acts on by permuting the 

tensor factors. Its action commutes with the action of GL(y). Therefore, for 
all TT G &ab we have that tt(x is a highest weight vector of type A. Indeed, 
Schur-Weyl duality (see e.g. |Pro071 Ch. 9, eq. (3.1.4)] or |(IW091 eq. (9.1)]) 
says that the set {tt^a | tt G Sab} is a generating set for H\N\/xi^^^^V). 

Let Tx denote the column-standard tableau of shape A, i.e., the standard 
tableau of shape A whose entries are consecutive in each column. For example, 
for A = (4, 2) we have 


1 

3 

5|6 

2 

4 



The group Gab acts on the set of tableaux with rectangular content (ab) x 1 
by sending the entries i to We can identify ttTa with tt^a- Indeed, the 

numbers in the columns of ttTx are exactly the indices of the position that are 
wedged with each other in ttCa- Since (x is a tensor product of special wedge 
products, the linear relations among all ttCx are given by the linear relations 
among the principal minors of generic matrices: The Grassmann and Pliicker 
relations. We can write these relations in terms of tableaux as follows: 

• T = —T' if T and T' differ in a single column by a single transposition. 
This is the Grassmann relation. 

• T = S, where the sum is over all tableaux S that arise from T by 
exchanging for some j and k the top k elements from the {j + l)th column 
with any selection of k elements in the jth column, preserving their vertical 
order, see |Ful971 p. 110]. These are known as the Pliicker relations. 

For example, from the Pliicker relations it follows that if T and T' arise from 
each other by exchanging two columns of the same length, then T = T'. We 
will use this in Section H31 

The Grassmann-Pliicker relations can be used to identify a set of tableaux 
that form a basis of HWVa({^“^ 14): Those are the standard tableaux, see |Ful97[ 
p. 110]. The algorithm that expresses a Young tableau with rectangular content 
(ab) X 1 as a linear combination over the set of standard tableaux is called 
straightening. 
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3.3.2 The inner and outer symmetrizations 

We now consider the canonical projection map s : ^ {^“Sym^y, which 

can be written as an element in the group algebra C[©Qb]: 


s := 


1 


E 

. 6 } 

0’2GS{b+l,6+2,...,2l)} 


(Ti(T2 • • • (Ja- 


f7a^&{b(a-l) + l,b(a-l) + 2,. 


We call s the inner symmetrization map. A tableau T is mapped to a sum of 
(6!)“ tableaux, as illustrated in the following example where a = 3 and b = 2. 


s 



2 

5|6 

3 

4 




2 

5 


2 

1 

5|6| , 

1 

2 

5 


2 

1 

5|6| 

3 

4 


3 

4 


4 

3 


4 

3 




1 

2 

6 

5L 

2 

1 

6151 , 

1 

2 

6 

5L 

2 

1 

6|5 


3 

4 


3 

4 


4 

3 


4 

3 



Here the pairs (1,2), (3,4), and (5,6) are symmetrized independently. 

The cases in which s{T) = 0 are those where positions are symmetrized that 
lie in the same column. In all other cases it can be shown (for example by tensor 
contraction) that s(T) ^ 0 and we use a short notation for the resulting sum of 
tableaux: In each summand, replace the first a entries I,..., a with 1, the second 
a entries a + 1,..., 2a with 2, and so on. Of course this gives the same tableau 
for each summand, so for example in the above example we obtain the tableau 
in m for each of the 8 summands. This tableau with rectangular content a x b 
shall now represent the sum s(T). If T is standard and s{T) yf 0, then s{T) is 
semistandard. The semistandard tableaux with rectangular content a x b form 
a basis of HWVA({^“Sym^I/). It is clear that on the space of tableaux with 
rectangular content a x b the Grassmann-Pliicker relations hold and that the 
straightening algorithm can be used to express a Young tableau with rectangular 
content a x b over the set of semistandard tableaux with rectangular content 
a X b. 

There is no such nice basis known for HWVA(Sym“Sym^T). This is an impor¬ 
tant open problem in representation theory: Find a combinatorial description 
of dim HWVA(Sym“Sym^I/), which is called the plethysm coefficient. We only 
know that 



T semistandard of type A with rectangular content a x b 


is a generating set of HWVA(Sym“Sym^lA). Since in our cases we have a < 6 
and hence |Ga| < 6! = 720, this is not a major computational problem. 

For a given Young tableau T with rectangular content a x 6 let us call 
g{T) := ^ {a^b)-symmetrized tableau. We will often write / = 

g{T). So the set of all (a, &)-symmetrized tableaux / of shape A is a generating 
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set of HWVA(Sym“Sym^F). An (a, 6)-symmetrized tableau shall be depicted by 
replacing the numbers by letters, thinking of the set of letters as an imordered 
set. So, for example. 



1 

3 

3 

0 

0 



A 

A 

C 

C 

B 

B 



3.4 Applying the map 

Using the tableau descriptions from 13.31 we can describe the map 

: HWVA(Sym“Sym''U) ^ HWVA(Sym'’Sym“U) 

in terms of tableaux. A generating set of the left hand side is given by the set 
of all (a, 6)-symmetrized tableaux of shape A, while a generating set of the right 
hand side is given by the set of all (6, a)-symmetrized tableaux of shape A. 

Proposition 11 (A combinatorial description of Given an {a,b)- 

symmetrized tableaux f of shape A, 4'^ ^(/) equals the sum over all tableaux 
S that are obtained by putting the numbers !,...,& on the positions of the A’s 
in T in any order, and putting the numbers 1,..., 6 on the positions of the B’s 
in T in any order, and so on. 

We highlight this combinatorial definition with an example. 




3,2 


A 

A 

C|C|\ _ 

1 

2 

1 |2| 1 

2 

1 

1 |2| 1 

1 

2 

2 |l| 1 

2 

1 

2 | 1 

B 

B 

) 

1 

2 


1 

2 


1 

2 


1 

2 





1 

2 

l|2| 1 

2 

1 

1 |2| 1 

1 

2 

2 |l| 1 

2 

1 

2 | 1 




2 

1 


2 

1 


2 

1 


2 

1 



= -4 


1 

1 

12 

2 

2 



where the last equality is the result of applying the straightening algorithm. 
Note that the map in Proposition [TT] is well-defined and that by grouping 
the summands correctly we can see that the image "^a,b{f) is a sum of {b,a)- 
symmetrized tableaux. 


Proof of Provostion \lll For an (a, 6)-symmetrized tableau / = X^o-g© 


T we 


take a summand T with rectangular content axb and interpret it according to its 
definition as a sum of tableaux with rectangular content [ab) x 1. Now we apply 
the map r, which means that we apply the corresponding permutation from &ab 
to the tableaux. We then turn summands into tableaux with rectangular content 
6 X a by applying the inner symmetrization map —?> {^^Sym“U. Upon 

close inspection one sees that we do not need to apply the outer symmetrization, 
as the resulting sum is already invariant under the action of ©&. □ 

Remark 12. In principle Propostion[TT]gives a combinatorial way of determin¬ 
ing the kernel of working with bases and using the straightening algorithm. 
The main obstacle is that the straightening algorithm is very slow for tableaux 
with many boxes. 
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3.5 Evaluating a function 

We fix p many random points Vi,...,Vp in Sym^E. To create a basis of 
HWV^(Sym“Sym^y) we choose random semistandard tableaux T of shape A 
and rectangular content a x 6, take its (a, &)-symmetrized tableau / = g(T) 
(which is an element of HWVA(Sym“Sym^y)) and evaluate / at the p points 
vi,..., Up. We choose more tableaux until we have found p linearly independent 
tableaux and discard any other tableaux. 

We proceed analogously for HWVA(Sym^Sym“F) with random points 
v[,...,Vp,, but here we store the points v{,... ,Vp, for later use (as described in 
Section 13.1|) . 

For this approach to work we need a method to evaluate a function given 
by an (a, 6)-symmetrized tableau at a point v € Sym^y. Analogously we need a 
method for the evaluation of (b, a)-symmetrized tableaux at a point v' G Sym“y. 
We will only describe the former one, the latter one being analogous. We will 
restrict our description to the case where u is a sum of a powers of linear forms, 
i.e., 

a 

v = J2^k)\ 

i=l 

where each U G V. We will always choose our points vj G Sym^y in a way that 
they satisfy this restriction. Moreover, we will choose li G 

Let T be a tableau with rectangular content axb and let / := g{T) . The prin¬ 
ciple of polarization and restitution explains the evaluation f{v) of homogeneous 
polynomials as a tensor contraction (where we interpret / G Sym“Sym^F C 

{g)“V): 

/W = (/W®“)- 

Since is symmetric under 6a, we see that 

/(u) = (T,u®“). 

Expanding we obtain 

a 

(T,u®“)= ^ (13) 

The following claim is the key insight of writing the evaluation in this way. 

Claim 14. Each summand in (nsi) is a product of determinants. 

Proof. Recall that since T has rectangular content a x 6 it is a sum of (&!)“ 
tableaux with rectangular content (ab) x 1. Note that the group (©&)“ stabilizes 
® ■ Therefore it is sufficient to contract 

where T is one of the {ab) x 1 summands of T, i.e., T = for some tt G &ab- 
The special form of (x (see equation (fTOll l implies that a contraction of irfx with 
a decomposable tensor is precisely a product of determinants, one determinant 
for each column of A. □ 
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4 Speeding up the algorithm 

The only problem with the algorithm stated in Section [3] is that it is too slow 
to do the required calculation in a reasonable amount time. In this section we 
present some crucial improvements that make the calculation feasible. 

4.1 The computation tree 

According to Proposition [TT] the naive way to compute 4^ ^ of an (a, b)- 
symmetrized tableau / is to loop over all ways to place the numbers 1 ,..., 6 
at the positions of the letters in / such that 1,... ,b gets placed on the A’s, 
1,..., 6 gets placed on the B’s, and so on. Many of the resulting tableaux will 
have a column in which two numbers coincide. By the Grassmann relation these 
correspond to the zero vector, so they do not contribute to the sum. We want 
to avoid enumerating these cases. A good solution is to place the numbers one 
after another, each time respecting the constraint that there should be no col¬ 
umn with two coinciding entries. This can be implemented with a computation 
tree and depth-first search. A priori this tree would have (6!)“ leafs, but this 
constraint prunes the tree considerably. 

A second important idea is to use the following common heuristic from com¬ 
puter science: While we are building the computation tree we are free to choose 
in which tableau cell we want to place the next number. We should always 
choose the cell where we have the fewest possibilities to place numbers. This 
heuristic tries to make the decision tree slim. 

4.2 Precaching 

If we do not want to compute 4^ {,(/), but are only interested in the evaluation 
'^ab(/)(^)! then at each computation leaf we compute a product of determi¬ 
nants. Before starting the algorithm we can quickly compute all the determi¬ 
nants that could possibly arise during the run of this algorithm and store them. 
This means that at each leaf of the computation tree we just perform a database 
lookup operation. 

4.3 Using symmetries 

If two columns ci and C 2 of / coincide, then the set of summands that we get 
when applying 4'o,b(/) can be partitioned into two sets: The ones where the 
top number in ci is smaller than the top number in C 2 and the ones where the 
top number in ci is larger than the top number in C 2 (equality cannot hold, 
since the two cells that we compare have the same letter in /). If we sum up 
either set of summands, by the Pliicker relations we get the same result. So we 
only have to look at one of the two sets. More generally, if several columns of / 
coincide, then we can enforce their top entries to be ordered. Moreover, for each 
top cell we can determine what its possible minimal and maximal value can be. 
For example, when calculating 46^6 and we have 4 coinciding columns Ci, C 2 , 
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C 3 , C 4 , then the top number of ci cannot be 4, 5, or 6 , as that would mean that 
the top number of C 4 would exceed 6 , which is impossible since we only place 
numbers 1,..., 6 . More generally, the (min,max) pairs for 4 coinciding columns 
are (1,3), (2,4), (3,5), (4,6) from left to right. This can readily be generalized 
to an arbitrary number of columns. These restrictions help to prune the search 
tree even earlier. 

These symmetry considerations give a large speedup factor if / has many 
coinciding columns. It turns out that if the plethysm coefficient is large, some 
of the tableaux have only very few coinciding columns. This is the main reason 
why we do not study kerb's,e completely. 

4.4 Parallelization 

The computation can be efficiently parallelized on several levels: First, each type 
A can be analyzed independently. Moreover, each tableau / can be evaluated 
at each point independently. Such an evaluation of a tableau at a point can be 
split up further by fixing a number of processors N and ensuring that in a fixed 
depth D of the computation tree the processor with number C G {0,..., — 1} 

only takes care of the Ith subtree if / = C mod N. 

4.5 Other speedups 

In principle the algorithm could be sped up further by choosing the random 
generic points a little less generic, so that some subsets oi {k \ 1 < i < p'} 
are linearly dependent. Then we could prune the tree even earlier. We did not 
implement this. 

It is also possible to not evaluate at a point, but to only contract with a 
symmetric tensor. This could speed up the computation if the number of rows 
is not maximal. We also did not implement this. 


5 Implementation and running time 

We implemented the above algorithm in C using the following external sources: 

• the BigDigits multiple-precision arithmetic library Version 2.4 originally 
written by David Ireland, copyright (c) 2001-13 D.I. Management Services 
Pty Limited, all rights reserved, 

• helpful C snippets from stackexchange for string handling, 

• code for ranking and unranking combinations from Derrick Stolee’s blog. 

We ran our code on the Texas A&M calclab cluster, which provides up to 
roughly 2000 CPUs for 9 hours each day, more hours on the weekend. 

While ker 45 5 can be determined on the cluster in a minute, the computation 
of the part of ker^g^e took several days and the computation for ker^g g ran 
for approximately 2 weeks. 
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